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The effect of discrete imperfections on the behavior of active transmission 
lines (i.e., lines with distributed gain) is considered. Two cases are studied: 

1. Lines with identical, equally spaced reflectors. The transmission and 
reflection gains versus frequency are studied as functions of the magnitude 
of the reflectors. Limits on the magnitude of the reflectors to guarantee sta- 
bility are investigated. 

2. Lines with rindom reflectors, having random position and /or magni- 
tude. The statistics of the transinission are studied; in particidar, the average 
value and the variance and covariance of the transmission are determined for 
small reflections. If the reflections become large enough, instability may oc- 
cur, and these calculations may become invalid. Stability of active distributed 
systems is studied in a companion paper} 

I. INTRODUCTION 

In the present paper we consider the theory of active transmission 
lines (i.e., lines with gain) with discrete imperfections. Both equally 
spaced, identical imperfections and random imperfections will be con- 
sidered. This study was suggested by R. Kompfner as a rough mathe- 
matical model for the effects of imperfections in certain types of optical 
maser amplifiers, in which the optical signal is reflected back and forth 
through the active medium on essentially nonoverlapping paths by an 
array of mirrors. A. G. Fox has suggested that this mathematical model 
will also provide a description of a one-dimensional active medium (e.g., 
maser) with (one-dimensional) random inhomogeneities. 

Consider an active transmission line that provides exponential gain to 
both forward and backward waves, and further provides distortionless 
amplification. The voltage (and current) then vary as 

e 2 — forward wave, 

+* K I . (1) 

e — backward wave, 
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Since the line lias gain, 



T = -a +j$. 



> 0. 



(2) 



(3) 



Since we assume distortionless transmission, the propagation constant 
/3 is related to the angular frequency u by 

/3 = u/v (4) 

where the velocity of propagation v is a constant independent of the 
frequency co. Further, the gain constant a is independent of a>. We 
may thus interpret /3 either as the propagation constant or as the 
normalized frequency. 

Consider a line with N discrete reflectors, as illustrated in Fig. 1. 
The wave traveling to the right at a distance z is denoted by W (z), 
the wave traveling to the left by Wi(z), as indicated in this figure. 
We take W (L k +) and W^L^) as the right- and left-traveling waves 
just to the right of the fcth reflector c k , W (L k -) and W x (L k -) as 
the right- and left-traveling waves just to the left of the fcth reflector. 

Each reflector is characterized by a scattering matrix relating inci- 
dent and reflected waves. Thus for the typical reflector illustrated in 
Fig. 2 we have 



"TTi(L*-)"| VW«{L k - 

.TF (L,+ )J " lWi(U+ 

s Js u Sl2 l 

[_Sn S22J 



(5) 



(0) 



r Ln — ""i 

r~ u 1 

t^LM^d pi pi 

c, c 2 c 3 ck-i ck c n ., c n 

Wo(o)-- 1 W (L k +)-^j Wo(Ln+)— *| 

! j ! 

|*-W,(<» ^ kw,(L k +) K— W,(L N +) 

1 Lk=2- 1 -L ' 

1 = 1 

Fig. 1 — Line with N discrete reflectors. 
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W (Lit) W (L k+ ) 
> 1 *- 

Ck 



W,(Lk-) W,(l_ k +) 

Fig. 2 — Typical reflector. 

The incident and reflected wave amplitudes are assumed normalized 
so that the power in any wave is simply the square of its absolute 
magnitude. For example, if the reflected wave is absent at the left of 
the obstacle in Fig. 2 the power in the incident wave is | W (L k — )| 2 ; 
similarly, if the incident wave is absent the power in the reflected wave 
is | Wi(L k — )| z . We make the following assumptions: 

1. The powers in the forward and backward waves are additive; 
for example, the total power P flowing in the -\-z direction at the left 
of Fig. 2 is given by 

P = \W (L,-)\ 2 - | TK,(X fc -)| 2 . (7) 

2. The reflectors are lossless, and consequently have unitary scat- 
tering matrices. For a reflector of a given magnitude there is a single 
arbitrary phase parameter in the scattering matrix; this phase has 
been chosen in such a way as to yield a scattering matrix for the obstacle 
of the following form: 



S = 



< lc| < 1. 






(8) 



c is a measure of the magnitude of the reflection; for c = the reflec- 
tion is zero and the guide is perfect, c is assumed to be independent of 
frequency, although this assumption is not compatible with physical 
realizability. We note that the matrix of (8) is correct only for co (or /3) 
> 0. For w (or /3) < the signs of the diagonal terms of the matrix 
must be changed, so that the various responses will be real, even though 
unrealizable; alternately, we may change the sign of c for negative co 
(or/3). 

Next consider the cascade connection of reflectors and ideal guide 
sections shown in Fig. 1. We require the wave matrix A corresponding 
to the scattering matrix of (8) for an obstacle. Referring to Fig. 2, 
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W (L k -) 






A k = 



VI 



i r i -jc k ~] 



(9) 



(10) 



The wave matrix for the fcth line section of length l k between reflectors 
Ca--i and en is given by 



W (L k -!+) 
|_Wi(I*-i+)J 



S'k 











'W (L,-Y 
W 1 (L,-) 



(11) 



Thus the matrix X k for the cascade connection of the fcth line section 
of length l k and the fcth reflector is given by 



[ 



Wo(Lt-i+) 

Wi(La-i+)J 



= X k 



Wo(L k +) 
lWi(L k +)l 

X,= 



r „+r/ t 



Vl - Ck 2 _+ i 7c*e 



-n* 






+n»- 



(12) 



The over-all wave matrix X for the line consisting of N sections in 
Fig. 1 is 



PTo(O) 
LTT,(0)J 

Setting 



= X 



W (L N +) 
\_Wi(Ls+) 



x = x,x 2 ■■ ■ x„ = Hx b . (13; 



[x u .rial 
.x'2i .122J 



(14) 



and referring to Fig. 1, the (complex) transmission and reflection 
losses L T and L R or corresponding (complex) gains G T and G R are given 
as follows: 



Lit — ■= — — 



L« = 



W (0) 
W (L„+) 

TTo(O) = *u 
Wi(0) Xn 



= .Tn 



(15) 

(16) 



W (0), Wi(0) and TF (Ljv+), the incident, reflected, and transmitted 
waves for the entire structure, are illustrated in Fig. 1. 

It has been necessary to state the above analysis in terms of wave 
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matrices that give the input as a function of the output (instead of 
vice versa) because the boundary conditions are known at the output. 
The output is assumed to be matched, so that in Fig. 1 

Wi(L„+) = 0. (17) 

In contrast, the reflection coefficient at the input is not known in ad- 
vance, and so it is not convenient to express the output -,- ; * ! 

\_yVi{L N -\-)j 

as a matrix product times the input tt// a \ • 

We consider below two cases of interest : 

(a) Identical, equally spaced reflectors, 

(b) Independent reflectors with random magnitude and/or position. 

II. IDENTICAL, EQUALLY SPACED REFLECTORS 

We now assume that all reflectors have identical magnitude and equal 
spacing. Setting 

c k = c, l k = I 
in (12), from (13) and (14) the over-all wave matrix becomes 



X = 



•I'll Xi3 
Xn £'22 



1 



r +ri +n-\N 

e —jc e 



By the usual methods we find : 



-rj 



(18) 



*» - (1 . ^ K+ _ jo) (K ^ N - *-"-*>' (19) 

Xn = (1 - #)Ak + - K.) (a+N ~ "-">• (20) 

a ± = cosh T/ ± Vsinh 2 Tl + c 2 , (21) 

+n -n 

K ± = -£±— -*=£-. (22) 

With the help of (15) and (16) the transmission and reflection gains or 
losses may be determined. 

Consider the various .r,, of (18), and in particular .r u and x 2 i of (19) 
and (20), to be functions of jftt, where we recall from (4) that is 
proportional to the angular frequency a>. We recall from the discussion 
following (8) that these results are valid only for positive frequencies, 
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> 0. The Xij have certain general properties of interest. First, we have : 
x ij \j(0l + »)] = (-D\P1, P ^ 0, (23) 

x*JU(r - fil)) = (-D^Vlffl, * S ». (24) 

Further, 

z.-,W#] = *</W ( 25 ) 

Equation (23) shows that a\, is periodic in the normalized frequency 0, 
of period 2ir/l. Equation (25) guarantees that the over-all response 
to a real input is real. Taken together, (23) and (24) show that the 
magnitudes of the losses | L r | and | L* | of (15) and (16) are periodic 
in of period r/l, and are symmetric about the points 01 = 0, x/2, 
t, 3t/2, • • • . Consequently in studying the magnitudes of these losses 
at real frequencies we need consider only the range ^ 01 ^ ir/2. 

Next, from (19)-(22) it might appear that the various functions 
Xij have branch points in the complex frequency plane because of the 
radicals in these equations. This is not true, however; a little study of 
these equations shows that the radicals really disappear for all (integral) 
N. Alternately, by considering the matrix multiplication of (18) it 
becomes clear that all the a:,-, are single-valued functions of r, and that 
no branch points can appear. 

We may thus determine the exact expression for the transmission 
or reflection gain via either (19)-(22) or direct matrix multiplica- 
tion in (18). However, we shall most often be interested in cases where 
the reflection parameter c is small in some suitable sense; application of 
perturbation theory to (19)-(22) greatly simplifies these relations 
and permits a useful interpretation of these results. 

Consider the radical in (21). If 

| c | « | sinh H | (26) 

then we may expand the radical in a power series and retain only the 
first correction term. Since 

| sinh IV | '"' = sinh 2 al + sin 2 01 ^ sinh 2 «/, (27) 

(26) will be satisfied for all if 

| c | « sinh al. (28) 

Therefore 

Vsinh' TJ + <? M sinh Tl + gj^jj . (29) 
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Then (21) and (22) become: 

**^*OTS- (30) 

K+ M -,' &r. E±H (3U) 

C 

*-*"*■ sera- < 31b > 

Substituting (30) and (31) into (19) and (20) and neglecting various 
small quantities, we obtain the following approximate results: 

*• = (I - c ?y,, «"" + «. (32a) 

'-Gara)^"""-"' (32b) 

jc -nsinhJVT/ /ooX 

»- (1 -^«. ' iShTT' (33) 

We make one further assumption, often used below, that the total gain 
in the absence of reflectors (c = 0) is large; i.e., referring to (2) and (3), 

e A ' a '»l. (34) 

Then (32b) becomes 

F ~ \2stihri) e ' e » L (35) 

So far we have ignored the question of stability; it is clear that such 
an active device can oscillate under some conditions. If the device does 
oscillate, our present results for loss (or gain) lack physical significance, 
for reasons discussed below. Instability can occur only if the gain func- 
tions of (15) and (16) have poles in the right-half complex frequency 
plane; if all poles of G r and G a are in the left-half plane the device 
will be stable. Since from (15-16) the poles of the G's are the zeros of 
Xu , we investigate the zeros of .r n as given by the approximate expres- 
sions of (32a) and (35). 

For c = 0, i.e., with reflections absent, the device will be stable, and 
consequently the zeros of .r„ lie in the left-half plane. It seems obvious 
on physical grounds that the device remains stable for small enough 
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values of \c\ , and will oscillate only when | c | exceeds some critical 
value. Assuming this to be true, we determine the conditions for stabil- 
ity by finding the minimum value of | c | for which a zero of x n appears 
on the real frequency axis, i.e., for some value of /3. 
From (32a) the zeros of x n occur when 

F = -1. (36) 
Equivalcntly, 

|F| = 1; (37a) 

Z.F = ±tt, ±3tt, ••• . (37b) 

Noting that 

sinh 2 17 = sinh 2 (- a + tf)l = (sinh 2 od + sin 2 fi) e~**, (38a) 

v = taT 1 £*# , (38b) 

tanh at 

where the principal value of tan -1 is implied, we have from (35)-(37) 
the following approximate relation for a zero of x n lying on the real 
frequency axis. 



•i 



/,' = c . e 2Nal e - j{2m - 2v) = -1. (39) 

4 (sinh 2 al + sin- 01) 

Thus 

Npl = <p + (tt/2) + imr; m = 0, ±1, ±2, • • • (40a) 



2 
C 2Nal 



4 (sinh 2 al + sin 2 01) 



- e 



= 1. (40b) 



^ is given by (38b). We now fix al and find the smallest value of | c \ 
for which (40) has a solution. Equation (40a), together with (38b), 
can be readily seen to have 2(N — 1) roots (/3/)> for < fil < 2w. 
For each of these roots there is a corresponding solution c = ± | c,- 1 
for (40b). It is obvious that the smallest of these | Cj \ corresponds to 
the smallest (#)> , which is that root lying closest to 01 = and which 

we denote (fil)i . 

For convenience we summarize the approximate results derived above 

in the present section. 
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*"" {!-*)»* e *" [1 +F] (41a) 



p _ ( C__ y -sat/ 

~ \2 sinh n/ e 



2ATaZ -j(2N0t 2*0 



4 (sinh 2 a£ + sin 2 0/) 
_i tan /3Z 

" = tan ssra- < 4lb > 

Conditions: 

I c | « Vsinh 2 al + sin 2 01 (41c) 

«* -I »l. (41d) 

The results of (41a) and (41b) will be valid for all if the condition of 
(41c) is replaced by the more restrictive condition of (42): 

| c | « sinh al. (42) 

The maximum value of the reflection coefficient magnitude | c | that 
yields a stable amplifier is given as follows, subject to the conditions of 
(41d) and (42) 

N(fl) t = tan"' jjjj^j + \ (principal value of tan" 1 ) (43a) 

I c | m „ = 2e~ Nal Vsinh 2 al + sin 2 (01) t . (43b) 

In deriving (43) we required that the results of (41a) and (41b) be 
valid for all 0. Consequently the more restrictive condition of (42) 
must hold; however, it is not obvious in advance that (42) will end up 
being satisfied in all cases. However, it is easy to show that this is in- 
deed so, so that the approximate limits on | c | imposed by the require- 
ment of stability are indeed given by (43), so long as (41d) is satisfied 
(i.e., the high-gain case). From (43a) we have 

(01)! < w/N. (44) 

From (41d) and (44) 

(0l)i « al (45) 

and consequently 

sin 2 (/3Z)i <K sinh 2 al. (46) 
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Equation (43b) thus guarantees that the more restrictive bound of 
(42) will always be satisfied in the high-gain case. 

The general behavior of the gain-vs-frequency (or 01) curve is readily 
seen from (41a) and (41b). In the second line of (41a) the first factor 
and <p vary slowly with 01, while the factor e~ j2m varies rapidly. The 
angle of F increases steadily as (31 increases from to 2tt; the magnitude 
of F is largest at fil = 0, tt, 2t, • • • , and decreases rapidly away from 
these points. Therefore the gain G r of (15) plotted vs/3Z (or frequency) 
will have an oscillatory behavior, with the magnitude of oscillation 
greatest near (31 = 0, tt, 2tt, • • • , and quite small elsewhere. The larger 
N, the more rapid will be the rate of oscillation. 

It is instructive to consider a few numerical examples. We consider 
the following two cases : 

20 log 10 e Nal s 20 log 10 e aLfl 

= 30 db, total gain in (i) and (it) below 
(i) 20 logio e al = 1 db, gain per section 

AT = 30, number of sections 
(180/tt) -(/3/)i = 4.05°, phase shift per section at oscillation 

| c | max = 0.00860, maximum value of reflection coefficient 
for stability 
(it) 20 logio e al =0.1 db, gain per section 
N = 300, number of sections 
(180/71-) • (flOi = 0.405°, phase shift per section at oscillation 
| c \ max = 0.000860, maximum value of reflection co- 
efficient for stability. 

The total gain in both cases is large, and hence | c | m . x has been com- 
puted by (43). The transmission gain G r plotted versus the normalized 
frequency pi for these two cases is shown in Figs. 3 and 4 respectively 
for several values of c. These results are computed by direct matrix 
multiplication [see (18)] rather than via (19)-(22) or via the approxi- 
mate results of (41). Figs. 3(a) and 4(a) show the gain vs normal- 
ized frequency for three values of [ c \ less than | c | ma x as well as for 
c = |c| ma x [computed via the approximate results of (43)], which 
corresponds to the limiting case of stability. It is readily seen how the 
device approaches instability as c approaches | c \ max . Figs. 3(b) and 
4(b) show computed curves of the "gain" versus frequency for a value 
of c greater than | c | max . Under these conditions the device is unstable, 
so that these curves have little direct physical significance; however, 
these curves do not look too different from the stable ones of Figs. 3(a) 
and 4(a). This should provide explicit warning against taking any such 
computed curve seriously without first investigating stability. 
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^f°- /3l, NORMALIZED FREQUENCY IN DEGREES 



Fig. 3 — Transmission gain vs normalized frequency for one-dimensional ac- 
tive medium with identical, equally-spaced reflectors. N = 30, number of sections; 
201ogi e a ' = 1 db, gain per section; total gain = 30db;c = magnitude of reflectors, 
parameter indicated on curves. 



A detailed picture of the behavior of these devices could be worked 
out in terms of the poles of the gain function in the complex plane. For 
small | c | the poles lie in the left-half plane. As \c\ is increased the 
poles move toward the ./-axis, causing greater oscillation in the gain- 
frequency curve. As \c\—*\c | mttx the closest pole touches the ./-axis, 
causing the gain to approach infinity at one frequency. Finally, as 
| c | becomes greater than | c | max this pole moves to the right-half plane 
and the "gain "-frequency curve becomes finite. As | c | increases further 
the first peak decreases, but the next pole approaches the ./-axis, so that 
the second peak increases, approaches infinity, and eventually decreases. 
The different peaks in the gain-frequency curve behave in a similar 
manner as the various poles cross the ./-axis in succession. 

Figs. 5 and 6 show similar curves for the reflection gain G« . G« 
approaches infinity for the same values of | c | and /3Z as does G r ; this 
must be so, since for the limiting case of stability, power must emerge 
from both ends of the device in the absence of any incident wave. As in 



272 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1964 

40,-* 



|W— 0.000711 



30 



O 50 

ui 
Q 



I- 




^r 



(a) 

STABLE 



— 


p.- ox 


00860 












(a) 

STABLE 


4--0.C 


00815 
















/ 




\ * 
















N 


V s 


s^^ 















fl Yt 


000 










(b) 

UNSTABLE 


M 


V v 


[/~\s* 



















I 


5 


* 


b < 




7 


J 9 



7p/3l, NORMALIZED FREQUENCY IN DEGREES 

Fig. 4 — Transmission gain vs normalized frequency for one-dimensional ac- 
tive medium with identical, equally-spaced reflectors. N = 300, number of sec- 
tions; 20 logio e al = 0.1 db, gain per section; total gain = 30 db; c = magnitude 
of reflectors, parameter indicated on curves. 

Figs. 3(b) and 4(b), the curves of Figs. 5(b) and 6(b) correspond to 
instability and hence lack direct physical significance. 

If the total gain in the absence of reflectors is not large, then the 
above results of (43) are not valid, and the approximate results of (41) 
are not valid over the entire range of permissible values of c. It is in- 
teresting to examine the exact computer solutions for one such case. 



(tit) 20 logio e al 
N 
20 logio e Na 

(180A).(/S7)i 

I ('- | max 



0.1 db, gain per section 

50, number of sections 

20 logio e aLN 

5 db, total gain 

5°, phase shift per section at oscillation 

0.065, maximum value of reflection coefficient for 

stability. 

Gain-frequency curves for several values of c are shown in Figs. 7 and 
8. The values of (/3/)i and | c | mB x given above have been determined 
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Fig. 5 — Reflection gain vs normalized frequency for one-dimensional active 
medium with identical, equally-spaced reflectors. N = 30, number of sections; 
20 logio &*' = 1 db, gain per section; total gain = 30 db; c = magnitude of reflec- 
tors, parameter indicated on curves. 



from these curves. As above, Figs. 7(a) and 8(a) show the transmission 
and reflection gains for the stable case, \c\ ^ \c | max , while Figs. 7(b) 
and 8(b) show the "gains" for an unstable case. The general comments 
given above for examples (t) and (it) apply also to this case. The 
approximation of (43), which was valid in examples (i) and (it) above, 
would have predicted (#)i = 3.37°, | c \ max = 0.0135 for the oscilla- 
tion conditions; this approximation is quite inaccurate in the present 
low-gain case, particularly for | c | max . 

.Straightforward calculation based on (18) or (19)-(22) in the peri- 
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Fig. 6 — Reflection gain vs normalized frequency for one-dimensional active 
medium with identical, equally-spaced reflectors. N = 300, number of sections; 
20 logio e al = 0.1 db, gain per section; total gain = 30 db; c = magnitude of re- 
flectors, parameter indicated on curve. 



odic case, or (12) and (13) in the general case, will of course always lead 
to some definite result for .r u as a function of frequency, whether or not 
the device is stable. However, only if we are assured that the device is 
stable will .r n have the desired physical significance of the steady-state 
loss function L r . If the device is unstable it will of course oscillate, and 
ultimately the linear behavior assumed here must break down. However, 
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Fig. 7 — ■ Transmission gain vs normalized frequency for one-dimensional ac- 
tive medium with identical, equally-spaced reflectors. N = 50, number of sections ; 
20 logio c al = 0.1 db, gain per section ; total gain = 5 db ; c = magnitude of reflectors, 
parameter indicated on curves. 

by demanding that the device be at rest at / = and examining the 
initial build-up of oscillation, the mathematical significance of x u may 
be examined in the unstable case. Suppose the device is initially at rest, 
and a sinusoidal input is applied at t = 0. The total response may be 
divided into a steady-state response, whose envelope is constant with 
time, and a transient response, whose envelope ultimately grows or 
decays exponentially with time in the unstable and stable cases re- 
spectively. The steady-state response is given by .r n in both cases. In 
the stable case, since the transients ultimately decay with time, only 
the steady-state response remains. In the unstable case the steady-state 
response retains the same mathematical meaning, but since the tran- 
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Fig. 8 — Reflection gain vs normalized frequency for one-dimensional active 
medium with identical, equally-spaced reflectors. N = 50, number of sections; 
20 logio e al = 0.1 db, gain per section; total gain = 5 db; c = magnitude of re- 
flectors, parameter indicated on curves. 

sient response grows exponentially with time, the steady-state response 
loses much of its physical significance. 



III. RANDOM REFLECTORS 

In the present section we consider active devices with reflectors having 
random position and/or magnitude; different reflectors are assumed 
statistically independent. Since the imperfections are random, the loss 
(or gain) is also a random variable, and we seek various statistics of the 
loss-frequency curve. The loss L T is determined from (12) -(15) ; we study 
the average loss and the second-order statistics of the fluctuations about 
the average, i.e., the variance and covariance of the loss fluctuations. The 
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form of (12)-(15) requires us to study the loss statistics rather than the 
gain statistics, which are of more direct interest. However, if the loss 
fluctuations about the average are small, then the loss and gain fluctua- 
tions will be almost identical (except for a change in sign), and their 
statistics will thus also be approximately identical. 

As discussed above, (12)-(15) yield the transmission loss L T only if the 
device is stable. If the device is unstable so that oscillation occurs, then 
the steady-state response L T given by (12)-(15) loses much of its physi- 
cal significance, as discussed in the previous section. The statistics of 
L T computed below are effectively averaged over all cases, so that these 
results will not be meaningful unless the probability of oscillation is so 
small that for practical purposes it may be ignored. Thus the results be- 
low are valid in the limit of very small reflections, in analogy to the per- 
turbation case of the previous section. In a companion paper 1 useful 
sufficient conditions guaranteeing stability are obtained; these stability 
conditions extend the range of validity of the present calculations to finite 
reflections. 

Three different statistical models of an active device with random 
reflectors are considered in the present paper: 
(i) random magnitude and spacing 

(ii) equal magnitude, random spacing 

(ill) random magnitude, equal spacing. 
Thus for case (i) in (12)-( In), c k and k will be random variables with 
appropriate distributions; we assume that the different r* and k- are 
independent random variables. In case (ii) the en are all equal to the 
same constant c , the k- are independent random variables. In case (Hi) 
the Ck are independent random variables, the lk equal to the same con- 
stant / . Case (ii) has been suggested by R. Kompfner as being appli- 
cable to certain optical maser amplifiers. 

In cases (i) and (Hi) we will assume that c* is symmetrically dis- 
tributed about 0, with a distribution narrow compared to 1. 

We assume in the present paper that h is always a large number of 
wavelengths, so that 

Ph » 2tt. (47) 

We further assume in cases (i) and (ii) that the distribution of l k about 
its mean is very narrow with respect to the mean, but wide compared 
to 27r//3. These assumptions are compatible with conditions existing in 
certain optical amplifiers to which these results might be applied. For 
certain calculations we need assume in addition only a smooth, sym- 
metrical distribution for l k about its mean. However, for certain other 
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calculations we must be more specific; here we will assume a Gaussian 
distribution for h , as follows: 

^'vt^-'^'- (48) 

where k is the expected value and o-, 2 the variance of h , 

k " «* (49) 

In accord with (47) and the discussion immediately following, we assume 
that 

2tt//3 « <tj « k ; cases (i) and (it). (50) 

Note that in case (m) 2* = k , as stated above, and o-, = 0. 

In the following work we make use of the Kronecker matrix product. 
For convenience we define this product and summarize some of its 

properties. 

Consider two matrices A and B with elements a,-> and by . The ma- 
trices A and B need not be square, have the same dimensions, or be 
conformable; their dimensions are completely arbitrary, so that the 
ordinary matrix products AB or BA may not exist. The Kronecker 
product, written as A X B, (as opposed to the ordinary matrix product, 
written as AB) is defined as follows: 



.1 X B = 



a,\\B anB duB 
chiB a-nB a^B 



(51) 



.4 X B has been written in (51) in partitioned form, with each sub- 
matrix consisting of a scalar element of A, a {i , multiplied by the entire 
matrix B. 

Kronecker products have the following useful properties: 

A X B X C = (A X B) X C = A X (B X C) (52) 

(A + B)X(C + D) -AXC + AXD + BXC + BXD (53) 

(,1 X B) (C X D) = (AC) X (BD). (54) 

As stated above, products without X's in (52) indicate ordinary matrix 
products, and the two matrices to be so multiplied must be conformable. 
Equation (54) may be extended to yield 

{At X Bi) (A 2 XB 2 ) ■■■ (A N X B N ) 

(55) 
= (AtA* ••• A N ) X (#i£ 2 ■■• B N ). 
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We now return to the results of Section I for the transmission of a 
general active device. From (13) we have (see Fig. 1) 



W'o(0) 
|Jf'i(0.)J 



= X&t 



x N 



TF (L A -+) 
lWi(L H +)Ji 



The output is assumed matched [see (17)], so that 

Wi(L*+) = 0. 

In computing the loss L T of (15) we might as well set 

W 9 (L N +) = 1, 

so that by (15) L T = TF (0); (56) then becomes 



L*i(0)J 



— X1X2 • • ■ Xn 



(56) 



(57) 



(58) 



(59) 



Now, in determining the average loss and the loss fluctuations about 
the average we are not particularly interested in the phase variations 
caused by the variation in total length, which may be large compared 
to the optical wavelength but is small compared to the average total 
length. Further, the variations in gain per section will also be small 
compared to the average gain per section. These considerations suggest 
the following transformations of (59), which remove these more or less 
irrelevant contributions to the loss and phase variations. From Fig. 1, 
the total length L N is 



L N = £ /, . 



(60) 



Next define £ T and (R as follows: 






L T = e +TL »'£ T , 


duT — 6 • LiJ- 


(61) 


WiiO) = e +TLN -(R, 


(H = e~ rz,Ar -PF,(0). 


(62) 



From (12) we define a new matrix Y k in terms of X k as follows: 

a a- — e • 1 k , 
where 



Y k = 



1 



vT^ 



or- 



1 



-jc k 

' e -2Tl k 



(63) 



(64) 
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Then from (60)-(64), (59) may be written 



m +TLif 



ft] 



,+ r 'i 



.+TI 



= e. 



+TL N 



Y 1 e +l,2 -Y 2 



•Y,Y 2 



A-Vls 



•'-a 



Cancelling out the e +TL " factor on both sides of (G5), 



((35) 



(66) 



where £ T is defined in (61), Y k in (64). 

Equation (66) is suitable for studying the statistics of the normalized 
loss £ T , which contains the essential information regarding the loss 
fluctuations of the device. The quantity <R has to do with the reflected 
wave at the input corresponding to a unit output wave, and will not be 
of further interest here. The factor e +TLN = e ~ aLN e iffLN removed from 
the unnormalized loss L T in (61) is of course a random variable, but 
for a given amplifier it has constant magnitude and delay. 

We now compute (£ T ), the expected value of the normalized loss £ T . 
Since the c* and l k are assumed independent random variables, the 
different Y k of (66) are independent random matrices in all three cases 
discussed above. Taking the expected value of both sides of (66), and 
noting that the different Y k have the same distribution, we have 



ra-Hj]- 



[67) 



where (Y) is obtained from (64) as 



(Y) = 



A 



Vi - c 2 / 



./ c _ 

- J wr= 



c : 



, ■ / c \ / -2ri 



I 1 \ . -2TI 



(68) 



Note that the independence of c k and h for a given k has been used in 
obtaining (68); the subscript A; has been omitted in the above relations, 
since the statistics of the different c k 's and of the different 4's are iden- 
tical. Finally, since we neglect the small variations in the gain per 
section, we may set 

(c- 2Vl ) « e ia %- m ) , (69) 
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where l u is given in (49) as the average length of the sections. Then (68) 
becomes 



(Y) = 



_-\ 
Vl - c 2 / 

c 



-J 



VT= 



I +>" e2al ° {vhi) v~> « 2a '° wh=i) <*"*"> 



(70) 



Now in cases (i) and (Hi) above we have 

c 



I c \ 



(71) 



since the distribution of c is assumed symmetric about 0. In cases (i) 
and (n) we have 



„-ffli 



(e~n « 0, 



(72) 



in view of the assumptions about the distribution of I. Consequently 
(70) becomes in the three cases: 



case (i) 
case (ii) 



(Y) = 



Wl - c 2 


\ 


1 


> 


1 


f" 1 — JCB 

_0 




Vl - Co 2 


> 


/ 1 

\Vl - c 


\ 

/ 


ri 

_0 e 2a 






(73) 



case (m) . 



From (G7) and (73) we have the following final results: 



/. 



-V 



(£t) = 



wr^T 2 / 



Vi-co 2 / ' 



cases (i) and (m) 
case (ii). 



(74) 



The result for case (ii) in (74) may be regarded simply as a special case 
of the results for cases (i) and (Hi). Since in cases (i) and (Hi) the 
distribution of c is assumed narrow compared to 1, we may in some 
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calculations make the following approximation in (74): 

1 



VT^ 



1 + h (c) , 



(75) 



where (c) is the mean square value of the magnitude of the reflection 
coefficient. 

Equation (74) shows that in all three cases the presence of random 
reflections has increased the expected value of the loss; further, the 
average loss is independent of |S and hence of frequency. Since (£ r ) ^ 0, 
if the deviations of £ T from its expected value are very small (as they 
must be in useful amplifiers), then we will have approximately 



(JB r > | « < | JBr | >• 



(76) 



This approximate relation permits us to estimate the variance of the 
magnitude of the loss, as discussed below. We note that 

|<£r>| ^ (\£t\). (77) 

Next consider the mean square value of the loss, ( | £ T | ) = (£r£r*). 
First note from (51) that 



[•] >< frl - 



£t£t 




rip i 2 1 


£ r (ft* 




£ r (R* 


(R£ T * 




(R£t* 


(R(R* 




j«r . 



(78) 



From (66), (55), and (78) we have 



\£ T \' 
£ T (R* 
(R£ T * 

l(R| 2 



= (Fi X Fi*)(F, X F 2 *) • • • (Y H X Y N *) 



(79) 



where Y k is given in (64). Taking the expected value of both sides of 
(79), again making use of the independence of the different Y k matrices 
and the fact that they have the same distribution, we have 



(80) 



where (Y X Y*) is obtained from (64) and (51) as shown in (81). 



~<|£r| 2 )l 




[l 


<c£ r (R*) 
<(R£r*> 


= (Y X Y*) N 






_<l<*l f >- 
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00 



*» -3- 3. 



I + 



^ '« r. 



+ 






1 + 



I 



•^-> '*■» 



K 

X 
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We again omit the subscript k in the above, since the statistics of the 
different c*'s and of the different fo's are assumed identical. 

We now apply the same assumptions used above to (81). As in (69), 
neglecting the small variations in gain per section leads to 

(e~ 2Tl ) « i*%~ m ) , 

{e- 2V ' 1 ) M e 2al "{e +m ) , 

4al\ ^ , Actio 



(82) 



(e )tt e 

where /„ as before is the average length of the sections [see (49)]. Further, 
we make use of (71) for cases («) and (in), and (72) for cases (i) 
and («). The resulting forms for (Y X Y*) differ in the three cases, 
but after some simplification the final quantity of interest, < | £7- 1 2 ) = 
(£ t £t*), is given by the following single relation in all three cases: 



(I £ T | 2 ) 

(I m| 2 > 

In case (n), we have in (83) 



/ 1 \ 



1 - c 2 



/ c 2 \ 



w 1 \ 



(83) 



/ 1 \ 



1 



/ 



\_ 



Co 



1 - Co 2 ' \l - C 2 / 1 ~ Co 2 ' 



(84) 



Equation (83) gives the desired result ( | £ r | 2 ) in terms of the nth 
power of a real matrix. The matrix power may of course be written out 
explicitly in the usual way, but for the sake of simplicity this will not 
be done here. Some numerical examples are worked out in the next 
section. The variance of the loss, denoted a £l 2 , is given by 



°zr - < I *r - <£*■> | 2 ) 

= (|£r| 2 >-|(£r>| 2 . 

The variance of the magnitude of the loss is given by 

-, £r , 2 = <[| £r I " < I £t I >] 2 > = ( I £t I 2 ) - < I £r I ) 2 

m < I £ T I 2 ) - I <£r> I 2 = <rzr\ 



(85) 



(80a) 



where the approximation of (86a) follows from (76). From (77) we have 

*, £ r I 2 $ **r % ■ ( 86b ) 

In these results ( | £ T | 2 ) is given by (83), (£ T > by (74); the approxima- 
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tion of (86a) should be good when t\ £t \/(£ t ) « 1. We see that for all 
thi-ee cases < | £ T | 2 ) and <t £t 2 are independent of /3 and hence of fre- 
quency. 

Finally we study the co variance of the loss £ T , denoted R £t (t), de- 
fined by 

R £t (t) = (£ T (P + r)£ T *(fi)) = R £t *(-t). (87) 

It will appear below that the expected value in (87) is indeed dependent 
only on t, and not on /3, within the approximations of the present treat- 
ment. If we regard the loss £ T (/3) as a random process, then the Fourier 
transform of R £t (t) yields the power spectrum of the random processes 
£ T (0)- R £t ( t ) thus gives information about both the dc and ac com- 
ponents of £ r (/3); of particular interest arc the mean square magnitude 
and the rate of fluctuation of the ac component of the loss. The total 
"power" (dc plus ac) P T of the random process £t(0) is 

Pr = R £r (0) = (\£r(fi) I')- (88) 

The dc "power" P dc of £ r (/S) is 

Pdc = R £ t(«>) = R£ T (-«>), (89) 

where the limits as r — > ± co exist. Both ac and dc "powers" arc neces- 
sarily pure real, and are of course independent of /3, since R £t {t) is 
independent of $ in general. Let us define the dc component of a given 
£ T {&) curve as 

1 r* 

£r dc = £ T (0) = lim 7 r J? £ T (H) d0, (90) 

.1/— oo iiB/L J—ht 

where the bar indicates an average over /3. Then it is easy to show that 
the dc power of (89) is also equal to 

Pdc - R £r ( « ) = R £r ( - ») = ( I £r d0 1 2 ), (91) 

where cC rdo is given by (90). Let us now define the ac component of a 
given £ T (P) curve by 

£tJ(3) = £r((J) - £ Uc . (92) 

Then the covariance R £t ( T ) OI the ac component £t- 11c (/3) and the ac 
"power" P.0 of the normalized loss £ T (0) are given as follows: 

R £t Jt) = (£ Tao (0 + r)jC r "*0S)) = R £r (r) - R £t («>) , (93a) 



Pac = {\£tJ&) f) = R £t (0) - R£t(°°) = Rstr (0). (93b) 
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For convenience we define the covariance of (R(/3) as an auxiliary 
quantity, although this quantity is not of present interest to us: 

RM = <<W + r)<R*(j8)>. (94) 

We have 



/fjB r (/S + r)l f 



£< 



(/3)1\ _ 


# £r (r) 
<JBr(/J+r)<R*(/5)> 
(<R(/3 + T)JBr*(/9)> 




'1 


= (Y(p + r) X F*GB)>* 










_0 



(95) 



(96) 



From (66), (55), and (95) 

ffj6r(0 
<JB r (/3 + r)(R*(/3)) 

<<R(/3 + r)JB r *(/3)) 

R<r(t) 

where we again make use of the independence of the different Y k and 
the fact that they have the same distribution. Using the various assump- 
tions given above in (69), (71), (72), and (82), and making appropriate 
simplifications in the different cases, we obtain the following final 
common result for cases (i), (ii), and (Hi): 



R £t (t) 



B«(r) 



/ 1 \ 



,. Io(r £ 



<•*■> 



1 - c 2 



(97) 



In addition to the usual approximations, we have used 

< e -^ +T) ') » (98) 

in cases (i) and (it) in obtaining the result of (97). This approximation 
implies that | r | « /3; i.e., we examine the covariance and hence the 
loss over only a relatively narrow (electrical) band. In the analysis we 
often use the quantity R £r ( « ), which gives the dc "power" [see (91)]; 
this is justified because the covariance computed from (97) will approach 
its asymptotic value R £t (<») for values of t satisfying the requirement 
| t | « )S. We assume the distribution of I is the Gaussian distribution of 
(48), and note that (e~ j2Tl ) is simply related to the corresponding charac- 
teristic function. Thus 



(99) 



f Note that this result justifies the approximations of (72) and (98) [subject to 
the condition of (50)]. A similar result for (e ri ), where r is complex, may be readily 
derived, and justifies the approximation of (69) and (82). 
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In case (Hi) we have a- t = in (99). Thus we have as our final result: 

t) 

(100) 



_«<r(t). 



/ 1 \ 

\l - c 2 / 



/ c 

\— 



'o /_ C \ „->2tJ „-2(t(t,)» 4a/ / 1 \ ->2tI -2(tct,)» 



\T= 



\5~=w 



/ 1 \ H 1 / g \ = Q> . 

\l - C 2 / 1 - Co 2 ' \l - C 2 / 1 - Co 2 ' 

o-i = 0; 



case (ii) 
case (ttt). 



(101) 



Certain general properties of R £t (t) are readily deduced from (100). 
First, R £t (t) is independent of /3 and dependent only on t, as assumed 
above in (87). Second, for t = 0, (100) becomes identical to (83), as 
it must. Finally, for r— > °o , we have incases (i) and (m) from (100) and 
(101) 



/ 1 



\r^" 



W-0 = 



1 - Co 2 / ' 



case (i) 

[also case (m) — see below] 

case (it). 



(102) 



^£r( °° ) i s rea l> as stated above. We recall from (91) that Rs, T ( °o ) is 
the dc "power" of £ T (P). The ac "power" is given by (93). 

Now, in case (Hi) the covariance Rj> t (t) is periodic, which implies 
that the random process £ T (0) is periodic; 4 however, this is obvious 
from the original formulation of the problem. Rz T ( °o ) no longer exists 
in the strict sense; the dc "power" is now the average value (over r) 
of R £t (t). It turns out that we may approach case (Hi) by considering 
case (t) and allowing a x to approach in (100). [This violates the con- 
dition imposed by (50) and used in the approximations of (72) and (98) 
and so the limiting process <t, — * is forbidden in some of the above 
results; careful examination shows that it is valid to allow cr t — * in 
(100).] Then R £t (t) does approach the limit of (102) as t-> co ; and 
so we take the first result of (102) as the dc "power" in case (ttt), as 
well as in case (i). 
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Iii general 



<£r(/3)>'- * Pac = (IMP) I), (103) 

tr^^P-m ( | £r B0 (/3) |"> . (104) 

However, in case (it) only — i.e., reflectors of identical magnitude and 
random spacing — (103) and (104) arc true with the 5^ replaced 
by =, as seen from (74) and (102). 

The matrix power of (100) is easily written explicitly in the usual way, 
but the results would be rather complicated. Numerical examples are 
worked out in the next section. 

IV. NUMERICAL EXAMPLE — RANDOM REFLECTORS 

Consider an optical amplifier with random reflectors of the type given 
in case (ii) of Section III: i.e., the reflectors have identical magnitude but 
random spacing. Assume: 

20 logioe"' — 1 db, nominal gain per section 

N = 30, number of sections 

20 logio e Nal ° = 30 db, nominal total gain. 

Fig. 9 shows the average normalized loss and the rms fluctuation of the 
normalized loss about its average value, plotted versus c , the magni- 
tude of the reflectors. As seen from example {%), Section II, instability 
is possible if | c 1 > 0.00860. Therefore the curves of Fig. 9 are solid 
for Co <0.00860, dotted for c > 0.00860. However, this is intended only 
as a symbolic reminder of the question of stability. We do not know 
whether or not instability can occur for | c 1 < 0.00860. Even though 
we know that instability can occur for | c J > 0.00860, the probability 
of instability might remain so small for some greater range of c that 
these curves would provide a useful approximation. In Ref. 1, Section 
VI, equations (122)-(131) we show that stability is guaranteed for 
I Co I < 0.00590, assuming that the maximum fractional variation in 
spacing of the reflectors [v in (124) of Ref. 1] is small compared to 1. 
This is indicated in Fig. 9. 

All of the above results have been independent of the precise distribu- 
tion of the Ik , the spacing between reflectors, except that the conditions 
of (47) and the following sentence must be satisfied. However, the 
co variance of the loss depends explicitly on the probability distribution 
of the /* . For our present example we therefore assume that the differ- 
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<£T(/3)>=-/p^=-/<|£ T (/3)| 2 )y 

AVERAGE NORMALIZED LOSS 




00086 



s 



0.002 0.004 0.006 0.008 0.010 0.012 0.0)4 0.016 0.018 0.020 



STABILITY GUARANTEED 
FOR VERY SMALL MAX- 
IMUM FRACTIONAL VARI- 
ATION IN SPACING 0.0059 
OF REFLECTORS 



C„, MAGNITUDE OF REFLECTORS 



Fig. 9 — Average normalized loss and rms fluctuation about the average for 
one-dimensional active medium with randomly spaced reflectors of identical mag- 
nitude. N = 30, number of sections; 20 logio V l » = 1 db, nominal gain per section; 
nominal total gain = 30 db. 

ent k are independent, with the Gaussian probability density given in 
(48)-(50). We further assume the following numerical values: 

(at/k) = 0.01, Co = 0.005. (105) 

Thus, the spacing between successive reflectors is accurate to about 1 
per cent, and the magnitude of the reflectors would guarantee stability 
in the equally spaced case of Section II. Of course a practical device 
would probably be built much more accurately, but the values in ( 105 ) 
are suitable for illustrating the general behavior. Fig. 10 shows the 
(complex) covariance #£,.„( t) of the ac component £r ac (/3) of the 
normalized loss for this case as a function of the normalized variable 
(Wtt)t, for < (/o/»)t < 4. Fig. 10(a) shows the magnitude 
| Bjj T(M) (t) I and Fig. 10(b) the phase Zfij»r M ( r ) + 58 ^ot; note that 
the linear component of phase has been removed in the plot of Fig. 
10(b). The covariance is seen to be approximately a damped periodic 
function of r; Fig. 11 shows a plot of the magnitude of the covari- 
ance at the points t = n(ir/h), which correspond closely to the max- 
ima of | R £t (r) | . 



290 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1964 
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0.0020 — 
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IT 




Fig. 10 — Covariance of ac component of normalized loss for one-dimensional 
active medium with randomly spaced reflectors. aJU = 0.01; c = 0.005, magni- 
tude of reflectors; N = 30, number of sections; 20 logio e al o = 1 db, nominal gain 
per section; nominal total gain = 30 db. 

We would expect some resemblance between the covariance of Figs. 
10 and 11, for reflectors with identical magnitude but random spacing, 
and the (nonrandom) case of Section II for reflectors with identical 
magnitude and spacing. For the nonrandom case we have seen that 
the loss is periodic; consequently the covariance will also be periodic, 
and will look something like that of Figs. 10 and 11 for the random case 
except that it will not be damped. Note that the large linear component 
-58 lor that has been removed from the phase curve of Fig. 10(b) 
implies that the power spectrum of the random process £ Tao (0) is con- 
centrated around the angular "frequency" -58 k ; this angular "fre- 
quency" corresponds to the rate of variation of the loss for two reflectors 
whose separation is equal to the nominal spacing of the two end re- 
flectors in the random case. 
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Fig. 11 — Approximate maxima of covariance of ac component of normalized 
loss for one-dimensional active medium with randomly spaced reflectors (see Fig. 
10). <T|//o = 0.01; Co = 0.005, magnitude of reflectors; N = 30, number of sections; 
20 logio e°'o = 1 dl), nominal gain per section; nominal total gain = 30 db. 



V. DISCUSSION 

The question of stability has been discussed for the periodic case at 
the end of Section II. There it is pointed out that these calculations are 
valid only if the device is stable, i.e., does not oscillate. The same is true 
in the random case. In the periodic case we can determine by calculation 
the limits of stability, and this has been done in the examples of Section 
II. Stability in the random case is studied in Ref. 1. 

Various higher-order transmission statistics may be calculated by 
methods similar to those used above, but the complexity of the calcu- 
lations increases with the older of the statistics. In addition, statistics of 
the real and imaginary parts of the normalized loss £ T may be readily 
determined by similar methods. 
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